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Discovery of supersymmetric (SUSY) particles at the Large Hadron Collider (LHC) has been 
studied for the models where squarks and gluino are much heavier than the lightest supersymetric 
particle (LSP). In this paper, we investigate the SUSY discovery in the models with degenerated 
mass spectrum up to jtilsp ~ 0.7mq. Such mass spectrum is predicted in certain parameter region of 
the mixed modulas anomaly mediation (MMAM) model. We find that the effective transverse mass 
of the signal for the degenerated parameters shows the distribution similar to that of the background. 
Experimental sensitivity of the SUSY particles at the LHC therefore depends on the uncertainty of 
" the background in this class of model. We also find that SUSY signal shows an interesting correlation 

f"*^ ■ between M e g and |? T which may be used to determine the signal region properly to enhance the 

' S/N ratio even if the sparticle masses are rather degenerated. The structure is universal for the 

£\j i models with new heavy colored particles decaying into visible particles and a stable neutral particle, 

dark matter. 
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I. INTRODUCTION: TRANSVERSE PHYSICS AND ITS NON-TRANSVERSE LIMIT 
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Minimal supersymmetric standard model (MSSM) is one of the most promising candidates of 



the physics beyond the standard model (SM) that may solve the hierarchy problem in the Higgs 



ysic 

sector The model predicts a set of superpartners (sparticles) whose charges are exactly the same 
to that of the SM partners but the spins are different by one half. 

These sparticles will be directly searched at the Large Hadron Collider (LHC) at CERN, which is 
a pp collider with 14 TeV center of mass energy. The LHC is scheduled to start in 2007. One of the 
o • interesting features of the MSSM model is the conserved R parity, and all sparticles are assigned to 
have odd R parity. Because of the conservation law, the lightest supersymmetric particle (LSP) is 
stable, and is assumed as a good candidate for the dark matter. At collider experiments, sparticles 
are produced in pairs, and a sparticle decay should produce at least one LSP. SUSY events at collider 
# ^ ■ experiments therefore have a distinctive missing transverse momentum. 

^ If supersymmetry (SUSY) is an exact symmetry, a particle and its superpartner have a same mass. 
Therefore SUSY must be broken spontaneously. The SUSY is considered to be broken in "hidden 
sector" , and the SUSY breaking will be mediated to our sector by mediation fields. The mechanism 
is already under severe constraints. For example, the mass of the sfermions with same charge must 
be universal so that they do not cause dangerous FCNC, or sfermion masses must be so heavy that 
correction from the SUSY sector is suppressed. This suggests that yet unknown symmetry /dynamics 
exists in the sector. The origin of the supersymmetry breaking and mediation mechanism may be 
understood indirectly by measuring masses of the SUSY particles. Therefore both the discovery and 
mass determination of the SUSY particle have been studied intensively in recent years 0, IE El • 

If SUSY is broken at some high energy scale and the boundary conditions are universal at the 
scale, strongly interacting (SI) SUSY particles are heavier than electroweakly interacting (EWI) 
SUSY particles. For example, gaugino masses follow the relation 

Mi : M 2 : M 3 ~ 0.4 : 0.8 : 2.4 (1) 

in the minimal supergravity model, therefore the gluino is six times heavier than the bino-like par- 
ticles. At the LHC, SI SUSY particles will be copiously produced in pp collisions, and subsequently 
decay into lighter EWI SUSY particles. 
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The gluino and squark decays are associated with jets with high transverse momentum (pt)- The 
transverse momentum is the order of the gluino and squark masses. Moreover, because the LSP is 
significantly lighter than the gluino, the LSP from the gluino decay also has a high p?. They would 
give a large missing transverse momentum to the SUSY events. In addition, decays of the EWI 
sparticles may produce high p? leptons. Events from the standard model (SM) processes do not 
have such high px particles. 

Motivated by these observations, following cuts are often applied to reduce the SM background 
events to the SUSY signal events |2j; 

• An event is required to have at least one jet with px > 100 GeV and three jets with p? > 50 GeV 
within \t}\ < 3, 

• The effective mass of the event must satisfy M e g > 400 GeV, where the effective mass is defined 
using the transverse missing energy ($ T ) an d the transverse momentum of four leading jets as: 

M eff = Yl PTi+^T (2) 
i=l,. .A 

If the event has hard isolated leptons, the effective mass may be defined as follows: 

M eS = PTi+ E Pti+$t, (3) 

i=l, ...4 leptons 

where sum of the lepton pt can be taken over the leptons with px > 20 GeV and \rj\ < 2.5 GeV. 

• The missing transverse energy must satisfy the relation: 

$ T > max(0.2M efr , lOOGeV). (4) 



• The transverse sphericity St must be greater than 0.2, where St is defined as 2A2/(Ai + A2), 
with Ai and A2 being the eigenvalues of the 2x2 sphericity tensor Sij = PkiPkj formed by 
summing over the transverse momentum of all calorimeter cells. 

To reduce the background further, hard, isolated lepton(s) may be required. These cuts are good 
enough to reduce the SM backgrounds from tt+n-jets and W(Z)+n-jets productions down to a 
manageable level, although the production cross section of the SM processes may be (9(10 4 ) higher 
than signal cross sections. While the SUSY production section reduces very quickly as sparticle 
masses increase beyond 1 TeV, but the signature peaks at higher M e g- where backgrounds can be 
ignored . Previous studies show that the squark and gluino with mass around 2.5 TeV can be found 
at the LHC in the minimal super gravity model (MSUGRA). 

In MSUGRA, the SM background after the cuts can be neglected safely. Then, the distribution of 
accepted events are also useful to determine the mass scale of SUSY particles. For example, the peak 
of the M e g distribution is sensitive to the squark and gluino masses. For the events with same flavor 
opposite sign dileptons, the invariant mass distributions, rriu, rriji, and rriju, are useful to reconstruct 
the SUSY particle masses m%o, m^o, m^o and mjo. 

Recently it is pointed out that a string inspired model based on the flux compactification (KKLT 
models) [5J] predicts a mass relation different from that of the MSUGRA @, 0, 0. This is called 
mixed modulus anomaly mediation (MMAM) model. This model has a volume modulas T and a 
compensator field of minimum supergravity model C as a messenger of the SUSY breaking. The 
SUSY mass spectrum depends on the ratio of the two SUSY breaking parameters Ft and Fc- The 
unification scale of the sparticle masses depends on the ratio. It is interesting that the unification 
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scale of the soft SUSY parameters can be much lower than the GUT scale in this model. They may 
unify even at the weak scale for a special choice of the model parameters. 

When sparticle masses are degenerated at the weak scale, we expect a reduced probability to have 
high pt jets, and smaller M e g and $ T for given squark and gluino masses. This means that the 
standard SUSY cuts reduce the signal events as well, and SUSY discovery is more affected by the 
SM background. 

Quantitative understanding of the SM background may be required in this case. Existing back- 
ground estimates have large theoretical uncertainty coming from the scale of the strong coupling. Re- 
cently, several groups has emphasized importance to include the matrix element correction j^, OH El 
in the previous parton shower estimate of the background, which significantly change the background 
distribution in the signal region. When overlap between the signal and background distribution is 
large, as it is expected in the degenerated SUSY spectrum, the uncertainty of the background must 
be taken seriously, and we also need to reconsider the cuts to reduce the background. 

The purpose of this paper is to illustrate the phenomenology of the degenerated supersymmetry. 
We take the MMAM model as an example. By changing the ratio Fq/Ft, the mass spectrum changes 
smoothly from MSUGRA-like one to anomaly-mediation-like one. In between, there are regions of 
parameters where squark, slepton, and gaugino masses are significantly degenerated compared with 
those expected in MSUGRA. The model therefore provides one dimensional parametrization from 
the "transverse signature" to its non-transverse limit. The investigation of our analysis may easily 
be extended to other SUSY/non-SUSY scenarios with ^ T signature, such as the universal extra 
dimension model or little Higgs model with T parity. 

This paper is organized as follows. In section 2 we describe the MMAM model and its mass 
spectrum with an emphasis on the region of parameter space where sparticle masses are degenerated. 
In section 3, we describe our Monte Carlo simulations. In section 4, we study how SUSY event 
distribution depends on sparticle mass degeneracy and compare it with the background. We find that 
the signal M e e distribution is quite similar to that of the background if m^o ^ 0.5rriq t g. However we 
find a certain relation among ^ T and M e g that is universal for the model with heavy colored particle 
that decays into the stable neutral particle LHC will miss the SUSY signature if m^o ^ O^m^g. 
Section 5 is devoted to discussions and conclusions. 



In this section, we briefly describe the MMAM model following the notation in p . In this model, 
all shape modulus and dilaton will be fixed by non-zero flux on a CY manifold in the Type IIB 
string theory. The low energy N = 1 Lagrangian of this model is given by unfixed volume modulus 
T, compensator field C, and gauge and matter fields W a and Qi as 



Here K e g = —3 ln(T + T*) + Zi(T + T*)Q*Qi, g c ^ v is the 4D metric in the super conformal frame. 
T = T + Ft9 2 and C = Cq + FcO 2 are the volume modulas and the chiral compensator superfield of 
N — 1 SUGRA, respectively. W e g = Wq(T) + \\jkQiQjQk is the superpotential for T and matter. 



II. THE MMAM MODEL 



A. Boundary conditions at GUT scale 




(5) 
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T dependent function and fa may be expressed as 

1 



(T + T*) 71 *"' f a (6) 

where rij is the modular weights and = 0(1) for matter fields located on D7 (D3) branes, and 
n—1/2 for matter living at brane intersections 0. 

In KKLT model, Wq = wq — Aexp~ aT , where the last term of Wq expresses the non-perturbative 
effect such as the gaugino condensation in D7 brane, which fix the volume modulus, Wq is the 
contribution of the flux. In addition to the N = 1 supersymmetric action, there are contributions 
from anti-D3 branes which break supersymmetry and uplift the potential from AdS vacuum to 
(nearly Minkowski) de Sitter vacuum. The term is expressed by a spurion operator depending on T 
and C, and minimum of the potential will be obtained by solving the effective N = 1 action and the 
lifting potential. 

The resulting theory is parametrized by Fq/Cq ~ m 3 / 2 and Ft/(T + T*) . The SUSY breaking 
terms are obtained by expanding the action by Ft and Fc- Here we define the soft terms as 

Lsoft = ~M a X a \ a - m^lQ,! 2 - AijkyijkQtQjQk + he (7) 
where y^ is a canonically normalized Yukawa coupling 



L 



soft 



fe-X'ZiZjZk 

They are explicitly written as functions of 7713/2 and R = m^j^iT + T*)/Ft as follows; 

M ° = (I + W) ^ 

2 (mi 1 jH, 1 <hi \ 2 
m * ~ [& + RdlnT 4d\nfi) m ^ 2 

Ajk = (s( m i + m j + ™>k) - -(li + lj + 7Jfe)) m 3/2 (9) 



R 



where m,- = 1 — m and 



^^sM 2 ?^-^ 2 ) (io) 

with C r = Y^ a Ta{r), namely for matter in the fundamental representation Cp = 4/3, Cp = 3/4, 
C = Y 2 . We only include the effect of the top Yukawa coupling y, therefore dq 3 = 1, d T = 2, 
dm = 3 and di = otherwise. 
The scale dependence of 7, is expressed as 

dfj, 327r 4ya 32vr 4 V2 y a J 



where 



1 Sign convention for the A parameter is such that the off-diagonal element of r mass matrix is — m T (A T + /jtan j3). 
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and 



dg b a 3 

= —<J • (13) 



d In fi 8n 

Here, D = £ r d r = 6, C 3 = £ r C r 3 = |, C 2 = |> C ' = §, and &' = 11, &i = 33/5 , 6 2 = 1, 6 3 = -3. 
Finally, 

(m, + m i + m fc )— y — 2^ ~T^9 a - ( 14 ) 



B. Mass spectrum in the MM AM model 



In Eq. ©, the highest term in 1/R is the contribution of modulas T, while the terms independent 
from 1/R is that of pure anomaly mediation. In 0, it is pointed out that the mass spectrum in the 
MMAM model shows a special feature if a ~ 1, where 



R 



a 



m(M p i/m 3 /2 



(15) 



This can be seen by investigating the low energy mass parameters as functions of a. For example, 
gaugino masses at the scale \i may be expressed as 



M a {n) 



m 3 / 2 
R 



1 - 



An 2 



Mi 



GUT 



{M Pl /m 3/2 ) a / 2 fiJ _ 



(16) 



in one loop level. When a ~ 2, the log term in the equation becomes at Msusy- Namely, if l a 
does not depend on the gauge group, the gaugino masses unify at the weak scale, rather than at the 
GUT scale. 

The mass spectrum in the matter sector depends on n, . If Yukawa couplings are non- vanishing 
only for the combination QiQjQk satisfying rii + rij + nk = 2, or if the effect of the Yukawa couplings 
in the RGE can be ignored, the scalar masses and trilinear couplings also unify at the same scale 
of the gaugino mass unification. This relation is satisfied for the choice n# = 1 and n matt er = 1/2 
and we call this choice of the boundary condition as model A. Even if this condition is not satisfied, 
squark and slepton masses tend to be close with each other at M SUSY when gaugino masses unify 
at low energy scale, because the gaugino loop corrections to the sfermion masses are roughly equal. 

In Table HJ we show the relation between R and a for a fixed gluino mass at the GUT scale of 
450 GeV. This corresponds to Mg(n) = a s (/i)/a(GUT) x M 3 (GUT) = 1.1 TeV. We can see that the 
unification of gaugino or squark and slepton masses occurs at R ~ 60. 



Q 


m 3/2 /7?(TeV) 


R 


0.26 x 10" 2 


0.45 


0.1 


0.30 


0.50 


10 


0.61 


0.56 


20 


0.92 


0.63 


30 


1.25 


0.73 


40 


1.58 


0.86 


50 


1.92 


1.06 


60 



TABLE I: The relation between a, 7? and m 3/2 /i?(TeV) for M 3 (GUT) = 450 GeV 
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For the following discussion, we calculate the low energy mass spectrum using ISAJET 2] version 
7.72 which solves the boundary condition given in section 2.1. The sparticle mass spectrum for the 
model A is listed in Table [HI As a increases, gaugino masses, and slepton and squark masses get 
closer, and the model shows the mass pattern different from that of MSUGRA. The numerical values 
roughly agree with those in 1 



R 


M 3 


Mi 


M 2 


A* 




m l 


me 


m 3/2 /R 


Q 


0.1 


1055 


184 


350 


700 


957 435 


354 


450 


0.26 x 10" 2 


30 


1045 


436 


536 


607 


913 


531 


476 


631 


0.92 


40 


1038 


573 


653 


545 


879 


578 


541 


729 


1.24 


45 


1034 


657 717 499 


852 


604 


576 


790 


1.41 


55 


1020 


882 


892 


339 


765 


671 


675 


951 


1.74 



TABLE II: Example of ISAJET solution of the low energy mass spectrum for model A, for different value of R. Here we fix 
the low energy gluino mass roughly constant by choosing M3 (GUT) =450 GeV and tan/3 = 10. All mass parameters are given 
in GeV. 



As can be seen in Table |Hj the higgsino mass parameter /z decreases as a increases. The [i is 
determined by solving the minimization condition of one loop higgs effective potential. In minimal 
supergravity model [R = limit), the higgs soft mass square at weak scale is driven to large negative 
value by stop/top loop, while the correction from gaugino loops is small because Mi, M 2 -C M 3 . The 
\x parameter is chosen to compensate the negative value to get correct gauge symmetry breaking, 
therefore /i can be as large as the stop mass in the minimal supergravity model. When a is large, 
M\,Mi 3> M3 at the GUT scale. The higgs soft mass parameters get large positive contribution 
from the gaugino masses. Top/stop loop effect is largely compensated by the gaugino corrections at 
the weak scale. As a result, M\ ~ fi for R ~ 37, and fi <C Mi, M2 for larger value of R. Although 
Mi and M2 get closer to M3 for R 3> 37, the mass splitting among SUSY particles increases again. 

We have seen in the model A that decrease of the \i parameter limits the mass degeneracy. In our 
study, we are interested in the model point where the mass splitting among the SUSY particles are 
the smallest. The \i parameter would be largest when the GUT scale value of the higgs soft SUSY 
breaking mass is the smallest. On the other hand, sfermion masses at GUT scale must be large to 
avoid the f\ or t\ LSP. We therefore also consider the model B, where uh u = riH d = 1 and n mat ter = 
as another example of the MMAM model. 

The mass spectrums of the models A and B for Ms(GUT) = 450 GeV and tan /3 = 10 are compared 
in Table IIHI The mass difference between squark (gluino) and LSP hits minimum when Mi ~ fi 
at R = R c . R c ~ 40 for the model A and 55 for the model B, where m^o/rriq = 0.55 and 0.70, 
respectively. As expected, the model set B has more degenerated mass spectrum at R c , because the 
fj, parameter is larger for the same gaugino masses. 

In the following sections, we discuss the discovery potential in the mass degenerated SUSY models. 
Phenomenologically important parameter for the hadron collider is the typical mass scale of the event. 
This can be expressed by the energy of the jet from q — > XiQ decay in the q rest frame, which is 
expressed as 

2p C M = (rn| - m|o )jm~ q (17) 

Because sparticles are pair produced, 2pcM is the typical order of the effective transverse mass M e g 
of the SUSY production process. We list the value of 2pcM for the model points in Table IIHI 



2 Note that ISAJET runs two loop RGE for all soft parameters, while the boundary conditions are calculated by the formula using one 
loop RGE. 
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set A 


set B 


R 






mu L {m- a )m % o 


2pCM 





995 (1055) 182 


961 


1041 (1061) 189 


1007 


10 


986 (1053) 246 


924 


1043 (1061) 248 


984 


20 


973 (1049) 326 


793 


1044 (1060) 330 


940 


30 


951 (1045) 426 


726 


1045 (1060) 434 


865 


40 


916 (1038) 507 


635 


1044 (1059) 569 


733 


50 


854 (1027) 426 


641 


1038 (1057) 713 


548 


55 


803 (1021) 335 


663 


1033 (1056) 721 


529 


60 


no EWSB 




1023 (1055) 700 


543 



TABLE III: The squark, gluino, and the lightest neutralino masses in model A and B. M3(GUT) = 450 GeV and tan/3 = 10. 



M eS oc 2p CM can be reduced by a factor of 2/3 to 1/2 for these models. In the following sections, 
we find that discovery of SUSY particles will be non-trivial in this region. It is worth noting that 
gaugino and higgsino are highly mixed when pcm is minimum. In this case, relatively large nucleon- 
LSP scattering cross section and smaller dark matter density are expected, possibly consistent with 
cosmological constraints. 

Another phenomenologically important aspect in terms of collider physics is the branching ratios of 
the gluino and squark into leptons. The leptons would be produced from the neutralino or chargino 
decays arising from the squark decays. For the model A, the decay q — > Xi (~ W, B) channels 
are open unless R > 50. The chargino and neutralino may decay into (s)leptons. Especially, if the 
decay channel q — > Xi ~^ I i} = e i f 1 ) ~ ¥ Xi is open, the large branching ratio into the golden mode 
q — > xqII is expected. For example, Br(% + — > fv) ~ 80% and Br(%2 ~~ ► ^ + ^ ) ~14% for the model 
A, M 3 (GUT) = 450 GeV, tan/3 = 10, and R = 37. For the model A, M 2 > mj if R > 20, namely, 
the decay channel is open especially in the degenerated region. For the model B, the squark and 
slepton masses are so heavy that the decay into slepton is always closed. 



III. MONTE CARLO SIMULATION AND RECONSTRUCTION 



Here we describe our event simulation method used in the next section. As explained earlier, 
SUSY mass spectrum is calculated by ISAJET [l2| which is interfaced to the HERWIG [13] event 
generator using ISAWIG Q. HERWIG generates hard processes, takes care of initial and final state 



radiations, and fragments partons into hadrons. 

To estimate event distributions to be measured at LHC detectors, we smear particle energies, 
identify isolated leptons, and reconstruct jets. We independently developed a fast detector simulation 
program 3 , which takes the following steps; 

1. Finding isolated leptons: If a lepton (e or fi) with Et > 10 GeV and \r)\ < 2.5 is found in an 



event record, we take a cone with a size AR = y (A77) 2 + (A0) 2 = 0.2 around the lepton. If the 
sum of Et of the particles in the cone (except the lepton at the center of the cone) is less than 
10 GeV, we regard it as an isolated lepton candidate. After the jet reconstruction described 
below, isolated lepton candidates with min(Ai?(/j)) > 0.4 are accepted as isolated leptons. 



3 Recently, a fine event simulator PGS (pretty good simulation) Il5l is also available. 
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2. Reconstruction of jets. We adopt a simple jet finding algorithm PYCELL, a subroutine in 
the PYTHIA package, with minor modification on its treatment of leptons and energy smearing. 
Namely, 

(a) We remove isolated lepton candidates for the jet reconstruction. 

(b) Particle energies are measured by the electromagnetic and hadronic calorimeters. We as- 
sume that the number of calorimeter cells in rj direction is 50 within \r)\ < 3 and 50 in 
direction, respectively. Transverse energy deposit Et in a cell is summed for electrons, 
photons, and hadrons as 

E T (EM) = E t+ E t 

electrons photons 

£ T (had) = E t, (18) 

hadrons 

which correspond to the hits in the electromagnetic and hadronic calorimeter cells, respec- 
tively. 

(c) Transverse energy deposit in each cell is smeared by Gaussian energy resolutions (Te t (EM) = 
0. 1 ^JE T (EM) [GeV] and cr£ T (had) = 0.5^J E T (had)[GeV]. After the energy smearing, we 
regard sum of the smeared Et(EM) and -Er(had) as the measured energy deposit in the 
cell. 

(d) We take the highest energy cell as the initiator and sum the energy deposits in the cell 
within AR < 0.4. We accept the cluster as a jet if J2 c \uster E T > 10 GeV. We repeat this 
procedure after removing the cells which have been already used for the jet reconstruction. 

The reconstruction process (1 and 2) is similar to the fast simulation codes used by the ATLAS and 
CMS groups. Note that we do not smear energy of isolated leptons and assume they are identified 
correctly, because the energy resolution and particle identification are excellent for leptons in both 
of the LHC detectors. We smear photon and non-isolated electrons energies with better resolution 
than that for hadrons, however we do not assume fine grained electromagnetic calorimeter in our 
simulation for simplicity. 

Our simulation does not take care of other detector effects, such as misidentification of leptons 
and hadrons, non-uniformity of detector responses (cracks), non-Gaussian smearing in the energy 
measurements. The simulation, however, reproduces the signal distributions of past simulation 
studies reasonably well. 

To check the validity of our independently developed simulation program, we compare published 
distributions of SUSY signals using the ATLFAST simulator [l6| with our simulation results. Figured] 
show distributions of the same flavor and opposite sig n dilepton invariant mass (m«) and the jet- 
lepton-antilepton invariant mass (m^u) at SPSla [17||, where ji is one of the two highest pr jets 
with m^ii < rrij 2 n. The distributions are important to reconstruct the cascade decay q — > x° — > 
I — > Xi- Contribution of background events with accidental two leptons can be subtracted using the 
distributions of e ± if ? events. The plots are produced by applying the same cuts as used in |4|; 

• M eS > 600 GeV. 

• at least 4 jets with p T1 > 150 GeV, px2 > 100 GeV, and Pt3,4 > 50 GeV. 

• two and only two leptons with pr(h) > 20 GeV and pr(h) > 10 GeV. 



• S T > 0.2 
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FIG. 1: The dilepton invariant mass distribution mu = m ee + m M(J — m eM (left) and rriju (right) distribution at the point SPSla 
by our fast simulation program. We produce 10 s SUSY events for the plot, corresponding to J dtC = 1.8 fb _1 . After the cuts 
and background subtractions, 445 dilepton events remain. 



• # T > max(0.2M eff ). 

We calculate the sphericity and $ T , based on the smeared energy deposits to the cells in \rj\ < 3 in 
addition to the momentums of isolated lepton candidates. 

We find that accepted number of events is consistent with the previous simulation studies. In 
Figure Q we show the m« and rriju distribution for 10 5 SUSY events (corresponding to J dtC = 
1.8 fb _1 ). The edge and the end point appear at correct positions. The uiu distribution has a 
sharper edge compared with that in pj, because we do not smear lepton energies in our simulation. 
We find that 445 events remain after the cuts and the background subtraction. In the report jij 
roughly 2.5 x 10 4 events remain after the cuts for / dtC = 100 fb" 1 . The acceptances of the two 
simulations therefore agree within la. 



IV. DISCOVERY OF THE DEGENERATE SUSY MODEL 
A. Mass degeneracy and signal distributions 

In this sub-section, we study the M e s distributions at the SUSY model points with degenerated 
mass spectrum. 

We first show the px distribution of the first jet in the left panel of Figure 121 Here we take the 
model B, M 3 (GUT) = 450 GeV, tan/? = 10 and vary degeneracy by changing R from to 55. As 
R increases, sparticle mass difference gets smaller, resulting in softer distribution. The peak 
positions are positively correlated with pqm which ranges from 504 GeV to 265 GeV when we change 
R from to 55. From the figure, we can read that the acceptance of SUSY events depends on R. 
For example, events with pt < 160 GeV is only 16% for the point with R = but it is 35% for the 
point with R = 55. This means that about one third of the events are rejected by the requirement 
on the pt of the first jet in the previous section for the latter point. 

The M e g distribution shows a similar behavior. In Figure 121 (right), we show the M e g distributions 
for the model B, with M 3 (GUT) = 450 GeV (solid histograms) and M 3 (GUT) = 650 GeV (dashed 
histograms), respectively. We apply the cut (A) and also require one hard isolated lepton with 
p T > 20 GeV and \r]\ < 2.5. 
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FIG. 2: Left) The p T distributions of the highest p T jet for model B, M 3 (GUT) = 450 GeV and tan ,3 = 10. 50000 events 
are generated for each points(i? = 0, 20, 30, 40, 45, 50, 55 from right to the left). Right) M e g distributions for 10fb _1 for model 
B with Af 3 (GUT) = 450 GeV R = 0, 30, 55 (solid histograms from right to left) and M 3 (GUT) = 650 GeV R = 0, 50(dashed 
histograms from right to left) respectively. 



Here we compare the distributions with R = 0.1 (MSUGRA like) and R = 30 to 55 (degenerated). 
It should be noted that, while the mass spectrums are considerably different, the power low of the 
distributions beyond their peaks are roughly the same. These high px events are originated from 
the collisions with its center of mass energy y/s much higher than the squark and gluino masses. 
Therefore the power low of the distributions only depends on the dominant luminosity function. 

The peak position of the M eS distribution has a direct correlation with the produced sparticle 
mass. In [3] it is found that the correlation between the SUSY mass scale defined as 

/ M 2 \ 

M^j SY = Msusy-v^ (19) 

where M x is the LSP mass and Msusy = (J2 o'i'm i )/{J2 The peak position of M eS is linear for 
MSUGRA and Gauge mediation model. We also find that the linear relation holds for the signal 
distribution. We do not provide the fit results here, because the existence of the standard model 
background is very important for this model, as will be discussed in the following subsections 4 . 

We now turn into the relation between E T and M e g for degenerated points. In Figure H3 we show the 
distribution for the model B with M 3 (GUT) = 650 GeV, tan/? = 10 and R = 0.1, 40, 50, respectively. 
The x- and y- axes are M eS and E T , respectively, in Figure El We see that the E T takes a significant 
part of M e g for degenerated points R = 40, and 50 if M e g- is higher, while, for R = 0.1 the E T for a 
fixed M e fj is broader in the plot. 

This distribution can be understood as follows. Suppose we have events with two uncorrelated jets 
with energy Ej et = Pcm- The dominant part of the cross section is squark/gluino pair production 
near the threshold for our case, so if they decay directly into two jets and two LSPs, the event 
kinematics are indeed of this type up to the boost to the beam direction. We now calculate missing 
energy E miss and effective energy M(eff) of the events. The momentum of the two jets and the 



4 The paper llSl also found the linear relation does not hold in general SUSY model. These are corresponding to the points where 
the dominant contribution to the total SUSY production cross section comes from lighter sparticles such as chargino, neutralino, and 
sleptons, which do not contribute to the 4 jet + missing Et signals. 
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FIG. 3: M e g vs $ T distribution for the model B, with M 3 (GUT) = 650 GeV, tan/3 = 10 and R = 0.1 (top left), 7? = 40 (top 
right) and R = 50 (bottom) respectively. 



missing momentum can be expressed as follows, 

Pi = (pcm,0,0,pcm) 
P2 = (pcM,Pcusine,0,pcucos9) 
Pmiss = (E miss ,-p CM sme,0,-p CM (l + cos9)) (20) 

where we take the direction of the p\ momentum as z-axis and events are in x-z plane. Then the 
missing energy E miss and effective mass M(eff) are defined as 

E miss = PcmV2 + 2cos#, 

M(eff) = 2 PCM + £miss, (21) 

where M(eff) ranges from 2pcM to 4pcM- 

The relation between E T and M e s is somewhat similar to that of -E m i SS and M(eff) as can be seen 
in FigureEl For each plot E T tends to be small when M e ff ~ 2pcM, where 2pcM = 1-36, 0.95, 0.7 TeV 
for R = 0.1,40,50, respectively. It increases linearly as a function of M e s up to M e s ~ 3pcM, and 
the number of events reduces quickly beyond M e g > 4pcM • The relation indeed should be true 
when two particle are pair produced and each of them decays into a visible particle and an invisible 
particle at hadron collider. 
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FIG. 4: The M e e distributions for non-lepton sample. The thick solid histogram is the background estimated in 0. Solid and 
dashed histograms are signal distributions for M 3 (GUT) = 450 GeV R = 0.1 (MSUGRA), and R = 37 (degenerate) in the left 
panel, and for M 3 (GUT) = 650 GeV R = 0.1 (MSUGRA) and R = 38 (degenerate) in the right panel, respectively. All plots 
are normalized for J dtC = 10 fb _1 . 



On the other hand, for the dominant background coming from W(Z)+ n-jets or ti +n-jets, the 
$ T /M e g should be smaller than those of the signal. This is because the SM events have many other jet 
activities which are not associated with the missing momentum productions. The background comes 
from highly boosted events, where ^ T is typically a fraction of the missing energy originated from 
high energy W, Z, and t. The canonical cut to reduce the SM background is ^ T > 0.2M ett -. If the 
background distribution in M e g- and plane is significantly different from the signal distribution, 
one can improve the S/N by applying an improved cut in the plane rather than the simple cuts 
M eS > 400 GeV and ^ T > 0.2M eS described in section 1. Because there are no background 
distribution study on the relation between M e s and ^ x , we do not attempt to improve the cuts in 
following sections and use the cut given in section 1. 



B. The zero lepton channel 

We now show the signal M eff distributions for the degenerated model. We start our discussion 
from the JvLr distributions of the SUSY events without high leptons. This channel has been 
studied in 0, H [H| with an emphasis on the model-independent reconstruction of the SUSY scale, 
as discussed in the previous subsection. 

In Figure EJ the signal distributions for the model A, with M3(GUT) = 450 GeV and 650 GeV 
are plotted for R ~ and R ~ 37. The signal distributions (solid and dashed histograms) are 
parallel to the background distribution (bars in the plots) beyond the peak of M e ff distribution. The 
difference appears only in the overall normalizations. As already discussed, this is because the signal 
distribution is determined mainly by the luminosity function, once M e g ^> 2Msusy- 

The background shown in the plots has a large uncertainty. The estimation has been made by the 
lowest order Monte Carlo generator ISAJET. The distribution is subjet to the scale uncertainty which 
is typically (9(30%). We also require at least 4 high p T jets, while the number of final state parton 
involved in the tree level diagram is much less than 4. The number of additional jets is estimated 
by parton shower (PS) approximation. The PS approximation is good to estimate distributions 
of jets collinear to the leading jets. Recently, several groups have emphasized the importance of 
the matrix element (ME) correction to the SM background processes. The matrix elements of the 
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FIG. 5: The signal M e g distributions of one lepton events. The solid and dashed histograms correspond to model A, R ~ 
(MSUGRA-like) and R ~ 37 (degenerated mass spectrum) respectively. M 3 (GUT) = 650, 750, 850 GeV from left to right. All 
plots are normalized for J dtL = 10fb _1 . 



diagram involvin g W 1 Z,t + n-partons are calculated by the generator devoted for the process with 
multi jets jjj UlL llflj] . Experiment al g roups are also working to take into account the corrections, 
for example, in the ATLAS group |l0| . the backgrounds from (it, W, Z)+n-jets are calculated up to 
n = (3,6,6) respectively using ALPGEN [11]. After the inclusion of the ME corrections, the SM 
background is increased by a factor of 2~4 for M e g > 2 TeV. Even in the highest bin in the plot, 
the expected background is above 1. This is caused by the increased high jets in the SM events. 
The improved background estimation is still preliminary and in the level of the lowest order, and 
subject of coupling scale uncertainty. 

We find that the signal rate is the same order of the background rate for the degenerated param- 
eters. If the uncertainty of the overall normalization of the background is 100%, no bound would 
be obtained from no-lepton channel by looking into M e g distribution only. However, as we noted 
previously, the structure of the signal distribution may be visible over the background distribution 
in the M eS and $ T plane. 



C. One lepton channel 

The significant fraction of SUSY production gives rise to events with isolated leptons (Pt > 20 GeV 
and \t]\ < 2.5). We show the signal distribution for the model A in FigureEl Here, we do not show the 
background distributions, because they are not given in previous literatures. In [lCj, the distribution 
of the SM background is given, which is roughly 1/20 of the background of the zero lepton channel 
with a similar shape. 

The histograms in Figure show signal distributions for M 3 (GUT) = 650 GeV, 750 GeV, 850 GeV 
from left to right, and R ~ for solid line and R = 38, 37, 37 for dashed lines, respectively. As can 
be seen in Table ITTT1 the minimum mass difference is smaller for the model B than that for the model 
A, because the higgsino mass \i is higher. Because of the degeneracy, the signal is more suppressed 
as can be seen in Figure El where we take M 3 (GUT) = 450 GeV (left) and 650 GeV (right), R = 55 
and 50 (dashed histograms), respectively. Overall normalization of the signal above M e g > 1.6 TeV 
reduces by a factor of 1/10 for for model B M3(GUT) = 650 GeV (Figure El right panel) compared 
to the model A (Figure 03 left panel) for most degenerate points. We conclude that one lepton signal 
is not accessible at the most degenerated points by looking into M e g distribution alone, assuming 
that SM background in the one lepton channel in [lQj . and a factor 2~ 3 uncertainty, even though 
the SUSY scale is relatively small for this point, rrig = 1.48 TeV The mass reach at the LHC will be 
revisited again in section 4.5. 
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FIG. 6: The M eff distributions for one lepton signal events for model B, M 3 (GUT) = 450 GeV (left panel), M 3 (GUT) = 650 GeV 
(right panel) and tan/3 = 10. Th solid histograms for R = 0.1 (MSUGRA like) and dashed histograms for R = 55 and 50 
respectively corresponding to degenerated mass spectrum. All plots are normalized for J Ldt — 10fb _1 . 




FIG. 7: The mu and vtiju distribution of signal for a degenerated SUSY point. 5 x 10 4 SUSY events are generated for model 
A, M 3 (GUT)= 450 GeV, tan/3 = 10 and R = 37. 



D. Two lepton channel 

For the model A, there are significant two lepton events in the \i ~ M\ region because I is 
lighter than the ly-like neutralino. We show the mjj and rriju distribution for the model A, with 
M 3 (GUT) = 450 GeV and R = 37 in Figure [7| Here we require that one of the two isolated leptons 
has pt > 20 GeV and the other has pt > 10 GeV. We also apply the standard cuts described 
in section I. The relevant mass spectrum for the most prominent endpoints are ma R = 901 GeV, 
mj fl = 869 GeV, m^o = 554 GeV, m^o = 499 GeV, and rrii R = 524 GeV. The predicted mu edge 

and rriju endpoint for the decay chain (Jr — > — > — y Xi is m u — 54.9 GeV and rriju = 308 GeV, 
respectively, consistent with Figure The edges and endpoints of other cascade decays are also 
visible in the plots. 

To select the jet from q decays, we have taken the one of two most highest px jets with smaller 
rrijii. If mass difference between rriq and m^o is smaller, we may found the smaller probability to 
find correct jets that arise from irtq decay We do not observe such effect in Figure [7| 
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The number of accepted events Ni (where i denotes number of isolated leptons with px > 20 GeV 
and \r]\ < 2.5) for this point are N = 11229, iVx = 5167, N 2 = 2645 for 5 x 10 4 generated events 
(corresponding to ~ 20 fb _1 ) . Number of two lepton events below the dominant edge is 2933 
after the background subtraction. Because M e g- distribution of the one lepton channel is not very 
prominent over the background, SUSY particles may be discovered in the two lepton channel for the 
degenerate region. 

Note that for the most degenerated SUSY spectrum, we have Mi ~ M 2 ~ fi so that all neutralinos 
are highly mixed. In such situation, dominant source of x° is QRi because qi has significant decay 
branching ratios_ into charginos. It is worth to point out that the relative weight of the decays 
Ql(r) — > X2 ~" > Il(R) ma y ^ e studied by l ooking into the charge asymmetry in rriji distributions, 



A = (rriji+ — rriji-) /(rriji+ + rriji-) |20j, |2l|. The charge asymmetry comes from the polarization of 
X2 arising from q, however as discussed in j22(], the effect tends to be suppressed if sparticle masses 
are degenerated. 

For the model B, the decay of \°i sleptons are closed for any value of R. For M 3 (GUT) = 
450 GeV and R = 40, we have N = 13433, iV a = 4412, N 2 = 642 for 5 x 10 4 generated events. 
While N and N± is about same as that of the model A, N 2 is reduced by a factor of 4. In addition 
no significant events with opposite sign same flavor leptons are seen after the subtraction of e ± n^ 
events. 



E. Discovery potential in one lepton channel 

In MSUGRA, the signal M e g distribution is clearly harder than the background. Especially, one 
can see the increase of the signal towards its peak, where the background is negligible. If the total 
M e ff distribution shows a bump or a clear change in the power behavior, one may claim discovery of 
SUSY particles without precise understanding of the background distribution. 

On the other hand, when SUSY mass spectrum is degenerated, the peak of the signal distribution 
shifts to a lower position where the background may be ~10 times higher. The signal distribution 
shows a similar power low behavior for M e R > M^ ak as can be seen in Figure 0] The bump structure 
may not be detected and precise understanding of the background would be required. Searching edges 
in mu distribution from the decay x\ ~ * ^ ~ * Xi may become more important, however, as we have 
discussed already, the signal rate is model dependent. We do not consider the possibility to improve 
S/N ratio by selecting signal region M e ff vs in this subsection. 

Table II VI shows the number of signal events after the cut for M e g intervals for 5 x 10 4 generated 
event (upper row) and number of events for 10 fb _1 (lower row) for both model A and model B. We 
take a moderate value of M 3 (GUT) = 650 GeV for the Table, which corresponds to irtg ~ 1.4 TeV. 

In the MSUGRA limit (R ~ 0), the background rate is negligible in the signal region. The number 
of events in M e R > 2.5 TeV is greater than 10 for R < 20, where the number of background events 
would be around 1 according to ^(j. However, the signal rate is reduced by more than a factor 
of 4 for R > 30 compared to that at R ~ 0. Another important question is if we have a flatter 
power spectrum in a certain M eS region. The number of background events is reduced steeply M eS 
bins in Figure EJ while we have nearly flat signal distribution between 1.2-2 A TeV for MSUGRA- 
like (R ~ 0) points. The ratio between the signal in the bins iV(1.6TeV < M e ^ < 2.0TeV) and 
iV(2.0TeV < M eff < 2.4TeV) ~ 3 : 1 for R = 30, 40 for the model A and for R = 30, 40, 50 for the 
model B, therefore the power law for the signal is similar to that of background in Figure 0]. We can 
describe the discovery potential in terms of the mass degeneracy rather than R. In the Table IIV1 
2p C M ranges from 1.4 ~ 0.7 TeV. (M 9 ~ = 1.4 TeV ). Once 2 Pc m reduces below 1 TeV, 70% of Mq, the 
separation between signal and background are not good in M e g distributions. Good understanding 
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Me ff (TeV) 0.8-1.2 


1.2-1.6 


1.6-2.0 2.0-2.4 2.4-3.2 


3.2-4.0 


a 


2p C M(TeV) 


model A 






R = 0.1 


448 


905 


1267 


999 


568 


69 


0.252pb 


1.35 




22.5 


45.6 


63.9 


50.3 


28.8 


3.4 






10 


752 


1201 


1378 


896 


483 


59 


0.245pb 


1.29 




36.8 


58.8 


67.4 


43.8 


23.9 


2.9 






20 


813 


1421 


1319 


676 


293 


17 


0.248pb 


1.20 




40.3 


70.5 


65.4 


33.7 


14.5 


0.8 






30 


1246 


1691 


1018 


336 


95 


9 


0.267pb 


1.04 




66.6 


90.4 


54.4 


18.0 


5.1 


0.5 






40 


2446 


2337 


758 


212 


63 


9 


0.314pb 


0.87 




152.3 


145.5 


47.2 


13.2 


4.0 


0.6 






50 


882 


960 


529 


198 


65 


4 


0.457pb 


0.90 




80.7 


87.8 


48.4 


18.1 


5.9 


0.4 






55 


173 


249 


175 


71 


26 


5 


0.666pb 


0.94 




23.0 


33.2 


23.2 


9.5 


3.5 


0.7 






model B 






R = 0.1 


450 


1108 


1405 


1192 


784 


96 


0.206pb 


1.42 




18.6 


45.7 


57.9 


49.2 


32.3 


4.0 






10 


649 


1267 


1585 


1145 


654 


67 


0.187pb 


1.38 




24.3 


47.4 


59.3 


42.8 


24.5 


2.5 






20 


836 


1528 


1645 


962 


432 


42 


0.174pb 


1.31 




29.1 


53.2 


57.2 


33.5 


15.0 


1.5 






30 


943 


1838 


1362 


561 


217 


26 


0.165pb 


1.2 




31.2 


60.7 


45.0 


18.5 


7.2 


0.9 






40 


1711 


1591 


711 


184 


86 


4 


0.158pb 


0.99 




54.0 


50.2 


22.4 


5.8 


2.7 


0.1 






50 


2428 


803 


215 


67 


25 


5 


0.154pb 


0.73 




74.9 


24.8 


6.6 


2.1 


0.8 


0.2 







TABLE IV: The number of events per M e ff bins in model A and model B. The upper rows show the events for 50000 events in 
our simulations, and lower rows show the event per 10 fb -1 . M3(GUT) = 650 GeV. 



of the background distribution and the detector effect is required to exclude the model in early stage 
of LHC runs / dtC = lOfb -1 . 

A factor 2 decrease of the signal is found for the model A between R = 50 and R = 55. In this 
model, squark mass decreases as R increases. While the production cross section increases as q gets 
lighter, the masses of heavier neutralino and chargino increase. They are heavier than q at R — 55, 
so that they would not be produced from q decay at all. The lighter neutralinos and charginos are 
still lighter than q but they are degenerated, so that only soft leptons are produced from the decay. 
The two lepton signal is also suppressed for this model point. This is an example that the rate of 
signal events with lepton is model dependent. 

The tendency of the signal reduction is same for the higher mass spectrum. The luminosity 
needed to discover the SUSY signal is larger at the degenerated SUSY points. For example, the 
production cross section is 43.6 fb for the model A with M 3 (GUT) = 850 GeV and R — 0, and 15.1 
events/10 fb -1 is expected in the region 2.4 TeV < M e s < 4 TeV, which may be enough to claim 
the discovery. On the other hand, for R = 37 only 3.8 events/10 fb -1 is expected, excluding the 
possibility of discovery with low luminosity. 
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V. DISCUSSION AND CONCLUSION 



In this paper, we study the discovery potential of minimal supersymmetric models with degenerated 
mass spectrum at the LHC. Such parameter regions have not been systematically studied in the past. 
SUSY studies have been mainly performed in the models where gluino and squark are significantly 
heavier than the LSP, i.e. in the MSUGRA model, gauge mediation model, and so on. 

We do not study the general MSSM model here, but take mixed modulas anomaly mediation 
(MMAM) models, which are parametrized by R oc Fq/Ft, tan/3, and M 3 (GUT). The parameter R 
is the ratio of the F terms of the volume modulas T in KKLT model and the compensator field C 
in N = 1 MSUGRA. This parameter R gives us one parameter description from the MSUGRA-like 
mass spectrum, where M3 3> Mi, to (moderately) degenerated mass spectrum. The models we have 
studied include the MSSM points where rriLSp/Mq^ ~ 0.7. 

When sparticle mass spectrum is degenerated, energy of a particle from squark and gluino decays 
is small relative to the mass scale. This is reflected in the distributions of the first px jet, $ T and 
M e ff. They tend to peak at smaller values than MSUGRA prediction for the same gluino and squark 
masses. The background from the SM processes increases rapidly for such low and low M e g 
region, therefore, the SUSY signal has to compete with large background. Indeed, the signal and 
background M e g- distributions are quite similar for degenerated parameters we have studied. We 
find the S/N ratio can be far below 1 for our most degenerated parameter points. Nevertheless, we 
found the signal distribution of degenerated SUSY mass spectrum show special universal pattern 
in M e fi and plane. This may help to determine the appropriate signal region and discriminate 
signal from background better. 

The discovery potential at the degenerated points therefore depends on our knowledge on the 
background distribution. The main background source to the SUSY productions is multijet final 
states involving W, Z and t. There are significant improvements in the lowest order estimation of the 
background recently; Matrix element corrections of the backgrounds are now included. Theoretical 
uncertainty comes from matching between parton shower approximation and matrix corrections, and 
higher order QCD corrections (scale uncertainty). 

Experimentally backgrounds are partly estimated by the experimental data itself. For example, 
overall normalization of dominant ti background can be determined by looking into the events with 
low where contamination of the SUSY events are small. For the degenerated mass spectrum, 
however, such calibration may also be affected. Note that the background estimation in the SUSY 
signal region require extrapolation using Monte Carlo, while there is a contamination in the region 
for the calibration in non-trivial manner for degenerated SUSY mass spectrum. We look into $ T 
distribution and find that the number of events with f£ T < 200 GeV is increased by factor of 2 for 
our most degenerated point (the model B with R = 55, Ma(GUT) = 450 GeV) from MSUGRA limit 
(# = 0.1). 

Even with the uncertainty discussed above, the discovery of supersymmetry in MSUGRA models 
would not be a problem up to ~ 2 TeV region. However, as we have discussed in this paper, this is not 
the case even for moderately degenerated SUSY parameters in mixed modulas anomaly mediation 
model. We stress the importance to include the background and its uncertainty to perform SUSY 
discovery and parameter studies for degenerated or general MSSM mass spectrum, while we find 
that they tend to be ignored in recent studies 23[ 24| . 
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